Optical Spectroscopy of Ha Filaments in Cool Core Clusters: 
Kinematics, Reddening and Sources of Ionization 

Michael McDonald^'^, Sylvain Veilleux^'^ and David S. N. Rupke^'^ 

ABSTRACT 

We have obtained deep, high spatial and spectral resolution, long-slit spectra 
of the H« nebulae in the cool cores of 9 galaxy clusters. This sample provides a 
wealth of information on the ionization state, kinematics, and reddening of the 
warm gas in the cool cores of galaxy clusters. We find evidence for only small 
amounts of reddening in the extended, line-emitting filaments, with the majority 
of filaments having E{B — V)< 0.2. We find, in agreement with previous works, 
that the optical emission in cool core clusters has elevated low-ionization line 
ratios. The combination of [O IIl]/H/3, [N II]/Hq;, [S II]/Hq;, and [O I]/Hq; allow 
us to rule out coUisional ionization by cosmic rays, thermal conduction, and 
photoionization by ICM X-rays and AGN as strong contributors to the ionization 
in the bulk of the optical line-emitting gas in both the nuclei and filaments. The 
data are adequately described by a composite model of slow shocks and star 
formation. This model is further supported by an observed correlation between 
the linewidths and low ionization line ratios which becomes stronger in systems 
with more modest star formation activity based on far ultraviolet observations. 
We find that the more extended, narrow filaments tend to have shallower velocity 
gradients and narrower linewidths than the compact filamentary complexes. We 
confirm that the widths of the emission lines decrease with radius, from FWHM ~ 
600 km s~^ in the nuclei to FWHM ~ 100 km s~^ in the most extended filaments. 
The variation of linewidth with radius is vastly different than what is measured 
from stellar absorption lines in a typical giant elliptical galaxy, suggesting that 
the velocity width of the warm gas may in fact be linked to ICM turbulence and, 
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thus, may provide a glimpse into the amount of turbulence in cool cores. In the 
central regions (r < 10 kpc) of several systems the warm gas shows kinematic 
signatures consistent with rotation, consistent with earlier work. We find that the 
kinematics of the most extended filaments in this sample are broadly consistent 
with both infall and outflow, and recommend further studies linking the warm 
gas kinematics to both radio and X-ray maps in order to further understand the 
observed kinematics. 

Subject headings: galaxies: cooling flows - galaxies: clusters - galaxies: groups - 
galaxies: elliptical and lenticular, cD - galaxies: active - ISM: jets and outflows 



1. Introduction 



The dense core of a galaxy cluster represents a unique environment, where the hot 
intracluster medium (ICM) is cooling most rapidly, feedback from the central active galactic 
nucleus (AGN) is most effective, and the brightest cluster galaxy (BCG) dominates the 
mass. This represents one of the few places in the Universe where large-scale cooling and 
feedback processes can be readily observed. Unlike isolated massive galaxies, where some of 
the energy injected into the interstellar medium (ISM) from the AGN often escapes into the 
low-density intergalactic medium (IGM), the denser ICM in cluster cores retains an imprint 
of this feedback in the form of bubbles (e.g.. 



Churazov et al. 


2001 


Fabian et al. 


2003 


). ^ 



of hot gas from the IGM onto massive galaxies is challenging to observe due to the very low 
densities, while such phenomena have been studied in depth for decades in galaxy clusters. 
The so-called "cooling flows" in galaxy clusters, which were once thought to be massive 
flows of cool gas on the order of 100-1000 M0 yr" 



see review by Fabian 1994), are now 

^ of cool 



understood to be considerably less massive, depositing on the order of 1-10 Mq yr 
gas onto the BCG (e.g. 



Voigt & Fabian 2004 Bregman et al. 2006). It is assumed that 



some form of feedback (e.g., AGN, gas sloshing, conduction) prevents the dense ICM from 
catastrophically cooling, allowing only a trickle of cool gas to accrete onto the BCG. 



There is a considerable amount of support for the "reduced cooling flow" model (Voigt 



& Fabian 2004) from observations at a variety of wavelengths. Cool core clusters (galaxy 
clusters for which the central ICM cooling time is much less than the age of the Universe) 
are well known to have an abundance of multi-phase gas in the central ~ 100 kpc which 
has been observed in O VI (e.g. 



2006), Ha (e.g. 



Bregman et al. 


2001 


Oegerle et al. 


2001 


Bregman et al. 


Johnstone et al. 


1987 




Heckman et al. 


1989 




Crawford 



eraL][T999| |Jaffe et al.|[2005j [Edwards et aL][2007l [Hatch et al.|[2007l [McDonald eFaL][20T0 
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2011a 


), and H2 (e.g 


., Edge 


2001 


Salome & Combes 


2003 


Jaffe et al. 2005 


Lim et al. 


2008 


Donahue et al.||2011 


). Additionally, evidence for star formation has been observed in the UV 


(e.g., 


McNamara & O'Connell 1 


989; Rafferty et al. 


2006 


McDonald & Veilleux ^ 


M09 


Hicks 


et al. 


2010 


McDonald et al 


201] 


-b), optical (e.g., Allen 1 


995; Cardiel et al. 


1995 


Crawford 


et al. 


1999 


Edwards et al. 


2007 


Bildfell et al. 2008 


, and mid IR (e.g., Hansen et al. 


2000 


Egami et al. 


2006 


O'Dea et al. 


2008 Quillen et al. 20( 


)8, hereafter MIR). In 


McDonald 


et al. 


(2011b 


), we showed that the star formation rates inferred from far-UV and mid-IR 



data suggest an efficiency of ~ 15% in converting the cooling ICM into stars, consistent 
with the global baryon fraction in stars. This star formation efficiency estimate, based on 
far-UV, Ha, and mid-IR data, is consistent with earlier work based on mid-IR star formation 



estimates by Egami et al. (2006) and O'Dea et al. (2008). This overwhelming evidence for 



cooling byproducts suggests that cooling flows are occurring, but on a reduced scale than 
was initially predicted. 



In McDonald et al. (2010, 2011a) we presented a sample of 33 galaxy groups and clus- 
ters with deep, high-spatial-resolution Ha imaging from the Maryland-Magellan Tunable 
Filter (MMTFQ [Veilleux et al. 2010). We reported on several new systems with extended, 
filamentary Ha emission, and produced higher-quality Ha maps for some previously known 
systems. By combining deep, archival Chandra X-ray Observatory (CXO) data with these 
new Ha images, we showed that the extended warm gas (Ha) is spatially coincident with the 
asymmetric, rapidly-cooling ICM for an ensemble of optically-emitting nebulae and that the 
Ha emission is always confined within the ICM cooling radius. This, taken together with a 
number of previously known correlations between the warm and hot phases (see review by 



Fabian 1994), showed that the Ha emission is intimately linked to the X-ray cooling fiow in 
both quantity and morphology. However, while we were able to determine the origin of this 
cool gas (cooling out of the hot ICM), we were unable to constrain the dominant ionization 
processes which produce the high Ha fiuxes. In an attempt to shed light on this remaining 
mystery, we have acquired deep, high-resolution spectra of multiple filaments and BCG nu- 



clei in a sample of cool core clusters drawn from our previous works (McDonald et al. 2010 
2011a|[b ). We present these data in §2, describing their collection, reduction, and analysis. 
In §3 we present the results of our analysis of the spectra, providing new information on 
the kinematics, ionization, and reddening of the warm, ionized gas filaments in cluster cores. 
In §4 we compare these results to various models in an attempt to ascertain the source of 
ionization in these filaments. Finally, we summarize our results in §5 and discuss future work 
which may shed further light on this topic. 



"'^ http://www.astro.umd.edu/ veilleux/mmtf/ 
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Throughout this paper, we assume the following cosmological values: Ho = 73 km s ^ 

MpC , ^matter 0.27, ^vacuum 0.73. 



2. Data Collection and Analysis 



Name 


z 


E(B-V) 


Source 


Slit PAs 


Ha Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Abell 0478 


0.0881 


0.517 


Keck - 


LRIS 


356.°4, 6.°4 


MIO 


Abell 0496 


0.0329 


0.132 


Keck - 


LRIS 


77.°4, 40.°4 


MIO 


Abell 0780 


0.0539 


0.042 


Magellan 


- IMACS 


110. °0, 72.°0 


MIO 


Abell 1644 


0.0475 


0.069 


Magellan 


- IMACS 


159.°0, 79.°0 


MIO 


Abell 1795 


0.0625 


0.013 


Magellan 


- IMACS 


350.°4, 344.°4 


MIO 


Abell 2052 


0.0345 


0.037 


Magellan 


- IMACS 


258.°1, 108.°0 


MIO 


Abell 2597 


0.0830 


0.030 


Magellan 


- IMACS 


14.°0 


Mllb 


NGC 4325 


0.0257 


0.023 


Magellan 


- IMACS 


33.°,5 163.°2 


Mlla 


Sersic 159-03 


0.0580 


0.011 


Magellan 


- IMACS 


86.°0, 161.°9 


MIO, 11a 



Table 1: Properties of our sample of 9 galaxy clusters with long-slit spectroscopy. The 
columns are: (1) Cluster name; (2) BCG redshift; (3) Galactic reddening toward cluster 



center from Schlegel et al. (1998); (4) Source of long-slit data; (5) Position angles of slits in 
units of degrees east of north; (6) Reference for MMTF Ha data. 

In order to investigate the various properties of the optically-emitting filaments in cool 
cores, we performed long-slit spectroscopy on a sample of 9 BCGs. This sample was chosen 
to include the most extended and luminous systems at Ha from our larger sample of 36 



clusters with deep Ha imaging from the Maryland-Magellan Tunable Filter (Veilleux et al. 
2010D as described in [McDonald et al.| ( |2010[ |2011a||b| (Table Q. Long-slit spectra were 
obtained over the span of 2 years at the Magellan and Keck telescopes, using the Inamori- 



Magellan Areal Camera & Spectrograph (IMACS; Dressier et al. 2011) and Low Resolution 



Imaging Spectrometer (LRIS; Oke et al. 1995), respectively. The setup we used allowed 
for a broad wavelength coverage, allowing full spectral coverage from ~5000-8000A with 
IMACS and ~3500-7500A with LRIS, and relatively high dispersions (IMACS: 0.385 A/pix 
@ 7000A, LRIS: 0.537 A/pix at @ 7000A), allowing for careful linewidth measurements down 
to the intrinsic instrumental linewidth (FWHM ~ 75 and 145 km s~^ for IMACS and LRIS 
@ 7000A, respectively). For each cluster, two slits were aligned along the filaments in an 
attempt to maximize the coverage of the optically-emitting gas. Additionally, one of the two 
slits was always forced to pass through the BCG nucleus. In the case of Abell 2597 only a 



- 5 - 



single slit position was acquired due to time constraints, and in NGC 4325 one of the two 
slits is slightly offset from the optimal position. The chosen slit positions are shown in Fig. 

m 

In total, 17 optical spectra of 9 cluster cores were obtained. These spectra were reduced 
using standard IRAF ( jhttp :7/iraf . netp packages for the reduction of long-slit spectra. 
This procedure involved removing the zero-point bias from and flat-fielding each exposure 
(CCDPROC), masking cosmic rays (LA Cosmic; http://www.astro.yale.edu/dokkuni7] 
[lacosmic/'), combining similar exposures (IMCOMBINE), removing spatial and spectral dis- 
tortions (TRANSFORM), removing sky lines (XVista; http : / /ganymede . nmsu . edu/holtz/J 
xvista/), and calibrating the spectra in both wavelength (IDENTIFY, FITCOORDS) and 
intensity (STANDARD, SENSFUNC, CALIBRATE). The wavelength calibration was de- 
rived from spectra of helium, neon, and argon arc lamps, while the photometric calibration 
was based on the spectrophotometric flux standards EG274, LTT3218, and G191-B2B. Spec- 
tra were visually inspected for contamination of the emission lines from the O2 atmospheric 
absorption features around 6870A and 7620A. 

From each spectrum, the redshift and velocity dispersion of the stars at the center of 
the BCG were measured using the Na D lines at 5890A and 5896A and the Mg I line at 
5178A wherever available. Next, the Ha and [N II] lines, which do not require deconvolution 
due to the high spectral resolution, were simultaneously fit using a combination of three 
Gaussians. It was assumed that all three lines had the same velocity dispersion and redshift, 
which yields 4 free parameters: v,ajj, Fna, ii]- In the few cases where multiple velocity 
components were visible, we fit two Gaussians to each line. The continuum was fit using a 
three-parameter model, which consisted of a linear component (2 parameters) and an Ha 
stellar absorption feature with the width and redshift fixed at the value measured from the 
Na D and Mg I lines. All spectra were visually inspected for contamination of the emission 
lines from the O2 atmospheric features at 6870A and 7620A. Of our 9 targets, only Abell 1644 
is redshifted such that the Ha+[N II] lines are affected by O2 absorption. For this system, 
we clipped the spectrum from 6870-6876A, which removed the strongly-absorbed part of the 
spectrum while leaving the majority of the emission lines intact. We note that, while the Ha 
line fluxes in this system are likely still contaminated and should be taken as lower limits, 
they do not appear as outliers in any of the observed trends and have H/3 fluxes and Ha 
FWHM consistent with the low observed flux, as we will show. 

Following the fits to the Ha+[N II] lines, the H/3, [O III], [O I], and [S II] lines were 
individually isolated and fit with Gaussians. For each of these fits the redshift and velocity 
dispersion were fixed to the values for Ha, while the line intensity and continuum zeropoint 
and slope were allowed to vary. Additionally, an H/3 stellar absorption feature was included 
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10 kpc 



Fig. 1. — Slit positions for all 9 cluster cores in our sample overlaid on Ha images from 
McDonald et al. (2010 2011a[ |b). The length of the slits have been artificially truncated for 
this plot to cover only optically-emitting regions. In all cases we have attempted to maximize 
the overlap between these slits and the Ha maps while forcing one slit to pass through the 
BCG nucleus. In the case of Abell 2597 we were only able to obtain one pointing due to 
time constraints. In NGC 4325 one of the slits is slightly offset from the ideal position due 
to a pointing error. 
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with a variable amount of absorption. Finally, a visual inspection of all 5 fitting regions was 
performed for each spectrum and spurious fits were re-run with more specific input parame- 
ters, allowing for a more reliable fit. All measured line intensities were further corrected for 



Galactic extinction following Cardelli et al. (1989) using reddening estimates from Schlegel 
eFaLl ((19981). 



3. Results 

These data, described in full in Appendix A, provide a wealth of information about the 
warm gas kinematics, ionization state, and amount of reddening in the cool cores of galaxy 
clusters. In this section, we first address the results of our spectroscopic analysis in the 
context of these three major topics, and then follow with a discussion of the implications that 
these results have on our current understanding of the generation of these optical filaments. 



3.1. Optical Line Ratios 



The use of deep, high spatial resolution, narrow-band Ha imaging to dictate the ori- 
entation of our long-slit spectroscopy allows us to characterize the properties of the optical 
emission in two spatial dimensions without the field-of-view limitations of 3-D spectroscopy. 
The broad wavelength coverage of these spectra allows the comparison of various optical 
emission lines, specifically the key diagnostic line ratios [N II] A6583/Ha, [O I] A6300/Ha, 
[S II] AA6716, 6731/Ha, and [O III] A5007/H/3. This combination of line ratios is tradi- 
tionally used to identify the source of ionization in emission-line regions ( [Baldwin et al. 



1981 Veilleux & Osterbrock 1987 Kewley et al. 2006). We show in Figure the first of 



these ratios, [N II] /Ha, which is often used to differentiate between star-forming regions 
and AGN, overlaid on Ha maps of individual cluster cores. Only 2/9 clusters (Abell 1644, 
Abell 2052) have [N II] /Ha strongly peaked on the nucleus, where we would expect AGN, 
if present, to dominate the ionization. Roughly 2/3 of systems have low (< 1.5) [N II]/Ha 
everywhere, while the remaining 1/3 of systems have high (> 1.5) [N II]/Ha everywhere. 



Interestingly, this is the same fraction of cool cores that appear to be star forming (McDon- 



ald et al. 2011b). For systems with both far-UV imaging from McDonald & Veilleux (2009) 



and McDonald et al. (2011b) and optical spectroscopy, we find a one-to-one correspondence 



between star-forming (Abell 1795, Abell 2597) and shock-heated (Abell 1644, Abell 2052) 
systems based on [N II] /Ha and far-UV/Ha flux ratios. The systems with the longest fil- 
aments (Abell 1795, Abell 2597, Sersic 159-03) have the lowest overall [N II] /Ha ratios, 
with the value of the ratio decreasing slightly at larger radius. However, these [N II] /Ha 




Fig. 2. — Spatial distribution of the [N II] A6583/Ha intensity ratio. The combination of 
narrow-band Ha imaging with long-sht spectroscopy allows us to create pseudo-2D spectra. 
In all panels, both the point size and color correspond to the magnitude of the line ratio, as 
described by the legend in the upper left. The cluster or central galaxy name is shown in the 
upper right, while the physical scale of the image is shown via a lOkpc scale bar in the lower 
right. The white cross represents the center of the optical (stellar continuum) emission. 
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Fig. 3. — Key diagnostic diagrams involving the reddening-insensitive line ratios [N II] 
A6583/Ha, [O I] A6300/Ha, [S II] AA6716, 6731/Ha, and [O III] A5007/H/3 for filaments 
(stars) and nuclei (circles) in each of our 9 systems. The grey points are galaxies from the 



Sloan Digital Sky Survey (Kewley et al. 2006) and show a clear separation into H II regions 



and AGN. The red line represents the extreme starburst limit and separates HII regions 



from AGN, while the blue line separates AGN into Seyfert and LINER classes (Kewley 



et al. 2006). Typical errorbars are shown in the bottom left corner of each panel. We note 



that, in general, the data for optical emission in cool core clusters tend towards elevated 
low-ionization lines. 
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ratios are significantly elevated above what one would expect for ongoing star formation 
(~ 0.3 for solar metallicity), but are roughly consistent with the expectations for starburst 
systems (~ 1.0; e.g., Kewley et al. 2006). Several prior studies have also found elevated 



low- ionization line ratios in the optical emission in cool core clusters (e.g., Voit & Donahue 



1997 Crawford et al. 1999). In order to understand these line ratios, we appeal to the full 



suite of line ratios first used by Veilleux & Osterbrock (1987). 



In Figure Isl we show the [O IIl]/H/3 line ratio as a function of [N II]/Ha, [O Ij/Ha, and 



[S II] /Ha. In these plots we show data from the Sloan Digital Sky Survey (SDSS; Abazajian 



et al. 2009) for emission-line nuclei, showing the separate regions occupied by star-forming 
galaxies, low-ionization nuclear emission-line regions (LINERs), and Seyfert galaxies. We 
find a systematic offset between the emission-line ratios in the warm filaments and the locus 
of points from the SDSS, towards lower [O IIl]/H/3 and/or higher [N II]/Ha, [S Ilj/Ha, and 
[O I] /Ha. We note that the grey points represent nuclei, while much of our data includes 
a stronger contribution from diffuse, ionized regions. The fact that excess blue light (e.g.. 



Crawford et al. 1999 Edwards et al. 2007 Bildfell et al. 2008), UV emission (e.g., Rafferty 



et al. 


2006 


Hicks et al. 


2010 


McDonald et al. 


2011b 


), mid- far IR emission 


et al. 


2008 


, and molecular gas (e.g.. 


Edge 


2001; Salome &; Combes 


2003) 



O'Dea 



is strong evidence for star formation, but Figure [3] seems to indicate that the situation is 
more complex. The [N II] /Ha and [S II] /Ha ratios in some filaments are consistent with 



the extreme starburst limit (Kewley et al. 2006), however the [O I]/Ha ratio is far too high 
to be pure star formation. Thus, it appears that pure star formation, whether ongoing or 
burst-like, is unable to explain all of the observed optical line ratios in the cool cores of 
galaxy clusters. In §4.1, we consider different combinations of ionization mechanisms which 
can explain the elevated low-ionization ratios (Figure |3| while still being consistent with 
the observations at other wavelengths that seem to strongly suggest the presence of star 
formation. 

The nuclei in our sample (colored circles in Figure |3]) have systematically low [O III]/H/3, 
consistent with LINERs. This classification is consistent with the relative weak luminosity 
of nuclei in cool core BCGs in both the optical and X-ray and with the narrow emission lines 



(Figure 10). However, only three nuclei (Abell 0496, Abell 1644, and Abell 2052) are clearly 
offset from the filaments in Figure [3j The remaining six nuclei are likely undergoing the 
same ionization processes as the filaments and, thus, may not be AGN, or may be composite 



objects (Kewley et al. 2006), despite their identification as LINERs in Figure 3 



The high spectral resolution of these data allows us to separate the [S II] doublet and 
determine the [S II] A6716/A6730 intensity ratio, which can then be used to probe the density 
of the line-emitting gas, as described in Osterbrock (1989). In Figure |4| the distribution of 



-li- 



the [S II] A6716/A6730 ratio is shown for both filaments and nuclei in our sample. The line- 
emitting gas in nuclei has typical densities between 100-300 cm~^, consistent with earlier 



work by Heckman et al. (1989). The more extended gas in the filaments appears to have 



overall lower densities, although we are unable to determine the average density due to the 
inability of the [S II] A6716/A6730 ratio to probe densities below ~100 cm~^. In general, the 
filaments tend to have N^. < 200 cm~^, potentially reaching much lower than ~100 cm~^ at 
large radius. 



Ng [cm-^] 

300 150 <100 




1.3 1.4 1.5 

[S II] A6716/A6730 



Fig. 4. — Distribution of the [S II] A6716/A6730 line ratio in the nuclei (upper) and filaments 
(lower) of 9 cool core clusters. The upper axis shows the inferred electron density of the 
optical line-emitting gas. At densities below 100 cm~^, the [S II] A6716/A6730 line ratio 
becomes ineffective and can only provide an upper limit. 



In Figure [5] we compare our measured [N II] /Ha line ratios to those from the literature 
in order to assess the quality of these measurements. There is generally good agreement 
between this work and the literature, suggesting that our measured ratios are reliable to 
within ~20%. Much of this scatter is likely due to different choices of aperture and slit 
orientation, and likely does not represent the absolute uncertainty in our measurements. In 
Figure M we compare our measured line ratios in nuclei and filaments to those measured in 



Perseus A (Hatch et al. 2006). In general, the line ratios measured at 98 different positions 



along the filaments in Perseus A agree with our measurements for the filaments in all 9 



-0.2 



-0.2 0.2 0.4 

log^o [N ll]/Ha„ 



iteroture 



Fig. 5. — Comparison of [N II] /Ha line ratios for this work to those taken from the hterature. 
Blue and red points represent extended and nuclear emission, respectively. Point types 



correspond to different references: Heckman et al. (1989) - filled circles, Melnick et al. (1997) 



- filled squares, Jaffe et al. (2005) - filled triangles. Hatch et al. (2007) - stars. The dashed 
line represents equality, while the dotted lines represent a scatter of ~20%. The significant 



outlier is Abell 2052, for which Heckman et al. (1989) had to decompose the blended [N II] 
and Ha lines. 



clusters. Furthermore, we find that the optical line ratios measured in Centaurus A (Farage 



et al. 2010) agree with both Perseus A and our sample. Overall, the filaments in cool cores 



tend to have elevated low- ionization ratios, with the [O I] /Ha ratio being very narrowly 
peaked at ~ 0.1. These line ratios are reminiscent of LINERs, but are measured far from 
the AGN, therefore suggesting a different origin. 



3.2. Gas Kinematics 

In Figure [7] we show Ha maps for each cluster core, with the line-of-sight velocity mea- 
surements overlaid. These velocities represent the line-of-sight velocity difference between 
the warm gas and the central (BCG) stellar component. We find that, in general, sys- 
tems with extended emission tend to have relatively smooth velocity fields (e.g., Abell 0496, 
Abell 1795, Abell 2597, Sersic 159-03), while compact systems and nuclei tend to have rota- 
tion signatures (e.g. Abell 0780, Abell 1644, Abell 0478). This observation is consistent with 
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-10 -10 -10 -10 
[N ll]/Ha [S ll]/Ha [0 l]/Ha [0 lll]/H/S 

Fig. 6. — Distribution of reddening-insensitive optical line ratios measured in Perseus A (top 
Hatch et al.||2006" ), nuclei (middle row) and filaments (bottom row) for the 9 clusters in 



row 



our sample. Note the overlap between the measured line ratios in Perseus A and in a typical 
filament. 



the work of Baum et al. ( 1990 ) and Baum et al. ( 1992 ) who found that emission-line nebulae 



in radio galaxies can be divided into three kinematic classes: rotators, calm nonrotators, 
and violent nonrotators. The systems which are classified as nonrotators (both calm and 
violent) tend to have relatively constant velocity fields and resemble the filaments in Figure 
[7[ while the rotators have strongly varying velocity fields, much like the nuclei and compact 
systems in Figure [7} This is more easily seen in Figure [8} which shows the line-of-sight 
velocity fields. In the central regions (< 10 kpc) the velocity fields in many systems show 
evidence for rotation, with a characteristic shape reminiscent of a rotating disk. At large 
radius (> 10 kpc), the velocity fields of the filaments appear to be relatively flat with typical 
line-of-sight velocities of ~300 km s~^ relative to the BCG. We do not find a correspondence 



between the kinematic class (rotator or non-rotator) as defined by Baum et al. (1990, 1992) 



and the ionization class (star- forming or shock-heated) described in §3.1 and McDonald et al. 



(2011b). We will discuss the implications of these kinematics further in §4.2. 
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Fig. 7. — Similar to figure [2| but showing the hne-of-sight velocity of the optical line- 
emitting gas relative to the brightest cluster galaxy. This figure shows that the more extended 
filaments (Abell 0496, Abell 1795, Abell 2597, Sersic 159-03) tend to have relatively smooth 
velocity fields with modest velocity gradients. 
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Fig. 8. — Line-of-sight velocity of the optically-emitting gas relative to that of the BCG 
for 9 clusters in our sample. Different slit positions are differentiated by point type, and 
adjacent spatial pixels are joined by solid lines. At small radii (r < lOkpc), there is evidence 
for rotation in several systems (Abell 0478, Abell 0780, Abell 1644, Abell 1795, Abell 2052) 



based on the measured v/a (Baum et al. 1992). Negative and positive radii correspond to 



positions along the slit with RA less than and greater than the BCG center, respectively. 
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Fig. 9. — Similar to Figure [2] but now showing the velocity width. In nearly all systems the 
velocity dispersion peaks at the position of the optical nucleus. In the more extended systems 
there is a trend towards decreasing velocity dispersion with increasing radius, suggesting that 
turbulence is larger in the cool cores. 
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Fig. 10. — Left: Similar to figure |8j but now showing optical linewidth as a function of 
radius for individual clusters. Right: Optical linewidth as a function of radius for all 9 
galaxy clusters. There is a strong correlation, such that the warm gas in the very center of 
the cluster is highly turbulent (FWHM ~ 400 km s~^) while the thin, extended filaments 
have very narrow lines (FWHM < 200 km s^^). This trend suggests that highly-elongated 
filaments can only survive in regions with minimal turbulence. Shown in blue is the stellar 
absorption line width in M87 from Murphy et al. (2011). 



In Figure [9] we show the velocity width of the filaments. In general, the velocity width 
peaks in the nucleus (FWHM ~ 400 km s~^) while, in the filaments, it is consistently lower 
(FWHM ~ 100-200 km s'^). The longest filaments (i.e., Abell 1795, Abell 2597, Sersic 159- 
03, NGC 4325) tend to have the narrowest emission lines. Curiously, Abell 2052, which 
harbors the most radio-luminous central AGN in our sample has the lowest central velocity 
width. Overall, the velocity width correlates well with the amount of turbulence one would 
infer from a visual inspection: long, thin filaments (e.g., Abell 1795, Abell 2597, Sersic 159- 
03) tend to have narrow emission lines, while short stubby filaments (e.g., Abell 0478) and 
regions with "disturbed" morphologies (e.g., Abell 0496, Abell 0780) tend to have broader 
profiles. 



Figure 10 shows the variation of the velocity width as a function of projected radius 
for each spatial element shown in Figure [9] There is an obvious correlation (Pearson R = 
—0.63) between the optical linewidth and radial extent of the filaments. This trend was also 



reported, albeit with lower significance, by Baum et al. (1990) and Baum (1992). This figure 
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shows compelling evidence that the most extended optical filaments are also those which 
are experiencing the least amount of turbulence. If we assume that the Ha-emitting gas 



is a product of the ICM cooling (e.g., McDonald et al. 2010), this may tell us something 



about the amount of turbulence in the ICM within the cooling radius. Beyond the inner 
lOkpc, the optical linewidth is < 200 km s^^, which is consistent with the observed upper 
limits on the ICM turbulence in cool cores from the XMM Reflection Grating Spectrometer 
(~ 200 km s"^; [Sanders et"ar||2011[ ). We note that the instrumental FWHM is ~ 75 km s~\ 



so our measurements are well above the lower limits. 

The fact that we see broader hues near the center of the cluster could also be due to 
the higher likelihood of chance superposition of filaments at smaller radii, where the number 
density of filaments is higher. While the kinematic data alone are unable to differentiate 
between a true increase in velocity dispersion and a projection of multiple narrow lines, we 
are able to offer evidence for the former by looking at the optical line ratios. In §3.1, we 
showed that the [N II] /Ha ratio is typically peaked in the nucleus, while in §4.1 we will show 
further that there is a correlation between the velocity width and the low ionization line 
ratios. These trends suggest that the broad line widths that are observed at small radii are 
real and are not strongly affected by projection. 



In Figure [TT] we compare our measured radial velocities and velocity dispersions to 
those from the literature. We find typical deviations of 1-2% and <30% in radial velocity 
and velocity dispersion, respectively. We note that the fractional error in relative velocities 
(i.e. velocity of gas with respect to the BCG) is much lower and is independent of our 
absolute velocity calibration. The larger deviations in the velocity dispersion measurements 
are most likely a result of differing apertures between this work and previous studies of the 
same sources. Overall, there is good agreement between our kinematic measurements and 
those from previous studies. 

These data tell an interesting story about the kinematics of the warm gas in the cool 
cores of galaxy clusters. In general, we find that longer filaments are less turbulent than 
shorter or morphologically peculiar filaments. The velocity width in the warm gas increases 
with decreasing radius, reaching a peak in the nucleus of ~ 400 km s~^, and a minimum 
of ~ 100 km s~^ beyond 10 kpc. The velocity fields are typically well-ordered along the 
filaments. We return to these results in 54.2 below. 
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Fig. 11. — Left: Comparison of BCG radial velocity measured for this work to those from 
the NASA Extragalactic Database for the same galaxies. This implies a typical error on our 
absolute velocity calibration of 1-2%. Right: Comparison of the velocity dispersion in the 
line-emitting gas for well-studied systems in the literature. Blue and red points represent 
extended and nuclear emission, respectively. Point types correspond to different references: 
Heckman et al. (1989) - filled circles, Melnick et al. (1997) - filled squares. Hatch et al. 



(2007) - stars. The dashed line represents equality, while the dotted lines represent a scatter 



of 35%. The large deviations are mostly dominated by the non-standard size of the apertures 
used in different studies. 

3.3. Intrinsic Reddening 



The amount of dust in the ICM is still uncertain. The abundance of relativistic particles 



should act to destroy dust grains on short timescales (Draine & Salpeter 1979). However, 
the fact that we observe star formation and, more importantly, mid-far IR emission (e.g.. 



O'Dea et al. 2008), suggests that some amount of dust may be shielded from destructive 



processes. Our data provide new estimates on the amount of intrinsic reddening in the optical 
filaments and nuclei in the BCGs of cool core clusters - something that is currently lacking 



in the literature for more than a few systems (e.g., Hu et al. 1985 Crawford et al. 1999). 
In order to determine the amount of intrinsic reddening, we began with the Hq;/H/3 line 
ratio, corrected for Galactic extinction and underlying stellar absorption (§2). We made the 
assumption that systems with Ha/H/9 > 2.85 (corresponding to case B recombination) had 
non-zero intrinsic reddening and determine the intrinsic E[B — V) assuming a dust screen 
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model and a standard extinction curve [R^ = 3.1; Cardelli et al. 1989). The reddening 
estimates here assume a foreground screen model and, thus, are underestimated by a factor 
of [exp(r) — l]/r if we assume that the line-emitting gas and dust are well-mixed. 



In Figure [12| we provide the 2-D distribution of reddening estimates based on {E{B — V)), 
corrected for Galactic extinction, for the 9 systems in our sample. Typical uncertainties in 
these estimates, assuming that the ~10-15% uncertainty in f/^^ dominates, is AE{B — V) ~ 
0.15. With the exception of a few outliers, the vast majority of systems appear to be 
mostly free of intrinsic reddening. The five most extended systems, Abell 1644, Abell 1795, 
Abell 2597, NGC 4325, and Sersic 159-03 all appear nearly reddening-free. There is a sig- 
nificant amount of reddening in the core of Abell 0780 (Hydra A), which is known to have 
a dust-lane at approximately the position of the reddening peak. Abell 0478, which has 
the highest Galactic E{B — V), seems to have a moderate amount of intrinsic reddening 
everywhere, but we cannot formally rule out a Galactic origin due to small fractional un- 
certainty on the large Galactic extinction value. Both Abell 0496 and Abell 2052 appear to 
have legitimate, non-zero reddening in the filaments. In the case of Abell 2052, the reddened 
filaments (north of the BCG nucleus) appear to be shock-heated (McDonald et al. 2011b) 
based on the high Ha/FUV ratio. 

We can also consider the distribution of reddening measurements in filaments and nuclei. 
Using the same regions as in Figure |3l we take the emission-weighted average E{B — V) in 



each region. The distribution of these values is shown in Figure 13 For the nuclei there is 
a relatively broad distribution in E{B — V) from 0.0-0.7, while in filaments the distribution 
peaks at E{B - V) = with a broad tail out to E{B - V) = 0.6. We find that ~70% of 
filaments have E{B — V)< 0.2, contrary to previous estimates based on randomly-oriented 
slits (e.g., Crawford et al. 1999) which found considerable amounts of optical reddening in 
cool core clusters. However, previous measurements contained only the nuclear emission 
and, thus, a comparison to our in-filament reddening estimates is unfair. We find good 



agreement between this work and Crawford et al. (1999) for the distribution of reddening 



in nuclei, as shown in Figure 13 The consistently small amount of reddening derived in 



the Ha filaments appears to be inconsistent with the evidence in the UV and mid-infrared 



Rafferty et al. 


2006 


O'Dea et al. 


2008 


Hicks et al. 


2010 



McDonald et al.||2011b ) , but is consistent with previous findings for the warm gas in Perseus 
A (E(B-V)~0.2; |Sabra et al.|[2000| and Centaurus A (E(B-V)=0.0-0.7; [Canning et al.|[20ll| ). 
We speculate that the harsh conditions surrounding the ionized filaments may prevent dust 



from surviving in all but the densest regions (Draine & Salpeter 1979). Finally, we note 
that these results are robust to our choice of Ha/H/3 = 2.85 as the dust-free limit and do 
not change dramatically if we instead choose Ha/H/? = 3.1, which refiects the conditions in 
shock-heated gas. 
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Fig. 12. — Similar to Figure |2|but now showing the intrinsic reddening, E(B-V), corrected for 
Galactic extinction. We note that many of the regions with locally high E(B-V) are coincident 
with low Ha surface brightness (e.g., Abell 0496, Abell 1795, Sersic 159-03) suggesting that 
they are more uncertain. 
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Fig. 13. — Distribution of optical reddening, E(B — V\ in the nuclei (upper) and filaments 
(lower) of 9 cool core clusters. The amount of reddening in nuclei is relatively flat over 
the range < E{B — V^< 0.4, while, in the filaments, the distribution peaks sharply at 
E{B - V") = 0, with a tail out to E[B - 1/) ~ 0.6. The distribution of E[B - V) in the 



nuclei matches well with that found by Crawford et al. (1999) (dotted line). The fact that 
the majority of filaments have little-to-no reddening suggests that dust does not survive long 
in the intracluster medium. The objects with E{B — V^) > 0.4 are Abell 0478, Abell 0496, 
and Abell 2052, as discussed in the text. These reddening estimates assume a foreground 
screen model and, thus, are underestimated by a factor of [exp(r) — l]/r if we assume that 
the line-emitting gas and dust are well-mixed. 



4. Discussion 

The data presented in §3 represent a significant step forward in combined spatial and 
spectral resolution and coverage over prior large spectroscopic surveys of emission-line neb- 
ulae in cool core clusters (e.g.. 



Heckman et al. 1989 Baum et al. 1990 Allen et al. 1992 



Crawford et al. 1999 Hatch et al. 2007; Edwards et al. 2009). Such a dataset allows us 
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to provide new constraints on the origin and ionization mechanism in these filaments. In 
the following subsections we discuss the implications of the results presented thus far and 
elaborate on their meaning. 



4.1. Sources of Ionization 



There are several possible sources of ionization in the cool cores of galaxy clusters. The 
most popular ideas have been: (a) the central AGN, (b) young stellar populations, (c) X- 
rays from the ICM, (d) heat conduction from the ICM to the cold filament, (e) shocks and 
turbulent mixing layers, and (f) collisional heating by cosmic rays. The lack of a strong 
gradient in [N II]/Ha with radius (Figure [2]) and the low [O III]/H/3 ratio in the filaments 
suggests that AGN do not contribute significantly to the ionization outside of the nucleus. 
The relative weakness of high-ionization lines (i.e., [O III]) suggests that ionization by ICM 
X-rays is also a small contributor. While scenarios (a) and (c) are relatively easy to rule out 
in the filaments, the remaining scenarios require a more quantitative approach. 



In Figure 14 we show the same data as in Figure |3| but now include model expectations 
for star formation, cooling plasma, shocks, conduction, and collisional ionization by cosmic 



rays. The model grids representing photoionization by young stars (upper panels; Kewley 



et al. 2001) show partial overlap with the data for filaments in the [N II]/Ha and [S II]/Hq; 
panels. However, these models fail to produce adequately high [O I] /Ha ratios, as noted by 
earlier studies (e.g., Crawford et al. 1999). The best match between the model and data is 



achieved if we use a model cloud with roughly solar metallicity and a low total ionization 
parameter (U). The location of this grid would move slightly if the IMF was altered, but no 
reasonable combination of ionization parameter, metallicity, and IMF would reproduce the 
observed [O I] /Ha ratios. 



In the second row of Figure 14, we show the expectation for condensing intracluster 



gas, originating at 10^ K (Voit et al. 1994). These models assume that the intracluster gas 
at T > 10^ K photoionizes the lO'^ K gas as it cools radiatively. In choosing the models to 
plot we make the assumption that the cool gas is ionization-bounded. This assumption is 
based on the fact that we observe strong [O I] and [S II] lines - species that are typically 



observed beyond the classical H II region (Hatch et al. 2006). These models reproduce well 



the observed [N II] /Ha and [O III] /Ha ratios observed in the filaments, but tend to slightly 
over-predict the [S II] /Ha and [O I] /Ha ratios. Interestingly, these models err in the opposite 
direction as the models for star formation, suggesting that a combination of the two may 
provide an adequate fit to the data. An important test for this model is whether the observed 
cooling rates can produce the high optical line luminosities that we observe. In [Donahue &: 




Fig. 14. — Similar to Figure |3} but with model predictions overlaid. From top to bottom: 
photoionization from young stars (lines of constant ionization parameter and metallicity; 



Kewley et al. 2001), self-ionization from a condensing plasma (lines of constant ionization 



parameter and metallicity; Voit et al. 1994), fast shocks (lines of constant speed and magnetic 
field strength; Allen et al. 2008 ), conduction (points correspond to different initial conditions; 
Boehringer & Fabian 1989), and coUisional ionization by cosmic rays ( [Ferland et al. 2009). 



In all panels, the point/color types are consistent with those in Figure [s] and the red dashed 
line is the extreme starburst limit of Kewley et al. (2006). 
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Fig. 14. — Continued. 



Voit (1991), the total Ha luminosity due to the condensing ISM is given as: 



L 



He 



-1 



Ho 



16 Mq yr 



-1 



(1) 



where Tj is the maximum temperature in units of 10^ K, fc is the fraction of the cooling 
radiation incident upon the photoionized clouds, and es is the fraction of incident radiation 
that reemerges in the Ha line normalized to 3%. We use 63 = 0.77, from Voit et al. (1994). 



In Figure 15 we show the Ha luminosity as a function of the X-ray-derived cooling rates 



for 14 clusters (red) from McDonald et al. (2010) and 7 groups (blue) from McDonald et al. 



(2011a). Even if we make the unrealistic assumption that 100% of the cooling luminosity is 
incident upon the cool clouds [fc = 1), the majority of systems are still much more luminous 
at Ha than predicted. Under these extreme conditions, the most luminous systems are 
consistent with still only < 20% of their Ha emission coming from the condensing ICM, in 



agreement with earlier results by Voit & Donahue (1990) using the same methods. Thus 



while the model of a self-radiating cooling flow is promising in terms of the predicted line 
ratios, the total fluxes are too low in the context of more recent, spectroscopically-derived 
X-ray cooling rates. However, the fact that we do measure non-zero cooling rates in the ICM 
means that self-radiation due to the cooling ICM is a contributing process and should not 
be ignored - a point we will return to later. 



In the third row of Figure 14 , we show model expectations for shocks of various speeds 



and magnetic field strengths from Allen et al. (2008). In general, the shock models cover 



the same range in [N II] /Ha, [S II] /Ha, and [O I] /Ha as the data, but slightly over-predict 
[O III]/H/3. The low-ionization ratios are best matched in the filaments by slow shocks 
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Fig. 15. — Total Ha luminosity versus X-ray-derived cooling rates from McDonald et al. 



(2010) and McDonald et al. (2011a). The red and blue points correspond to galaxy clusters 



and groups, respectively. The dashed lines show the predicted Ha luminosity assuming that 
the flux is due to self-ionization by a condensing hot plasma, for covering factors of 25% 
and 100%. The dotted line is five times the expected Ha luminosity for a covering factor 
of 100%. For the most extended and luminous filaments, this process contributes, at most, 
20% of the Ha flux, assuming an unrealistically high covering factor. 



(f ~ 100 — 400 km s~^), and in the nuclei by fast shocks [v ~ 800 km s^^). This makes 
qualitative sense, since the line widths are considerably broader in the nucleus than in the 
filaments. To further investigate the role of shocks, we consider the velocity dispersion as 



a function of the [N II]/Ha, [S II]/Ha, and [O I]/Ha ratios in Figure 16, We find that the 
systems with the weakest low-ionization line ratios have velocity dispersions consistent with 
those found for luminous and ultraluminous infrared galaxies (hereafter LIRGs and ULIRGs, 



respectively), which are likely experiencing a mix of shocks and star formation ( Veilleux et al. 



1995 1999 Monreal-Ibero et al. 2010). The systems which have low ionization line ratios 



considerably larger than those found in ULIRGs and LIRGs were identified earlier as non star 
forming based on the optical line ratios and far-UV emission. In the right panels of Figure 
16, we show that these non star forming systems have [N Il]/Ha and [S II]/Ha ratios which 
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are more strongly dependent on the velocity dispersion than their star forming counterparts. 
One possible interpretation of this result is that the filaments and nuclei are ionized by a 



mix of shocks and star formation. In this scenario, the two lines in Figure 16 represent two 
extremes: the shock-dominated case (steep line) and star formation-dominated case (shallow 
line). We will further investigate this scenario in the next section. 



In the fourth row of Figure 14 we show the expectation for heat conduction along the 



boundary between the hot ICM and the cool filaments from Boehringer & Fabian (1989). 
This mechanism produces either too high [O III] /Ha ratio or too low [N II], [S II], and [O I] 
ratios, depending on the choice of conditions. While the presence of soft X-ray emission 
coincident with optical filaments is evidence that conduction may be occurring (e.g., Nipoti 



& Binney 2004 ) , it cannot be the dominant source of ionization of the optical filaments on 



the basis of the observed optical line ratios. 



Finally, the bottom panels of Figure 14 shows the predicted line ratios for collisional 



ionization by cosmic rays (Ferland et al. 2009). While this model accurately predicts the 



[S II] /Ha and [O I] /Ha ratios, it underpredicts the [N II] /Ha and [O III]/H/3 ratios by 1 
and 3 orders of magnitude, respectively. In the absence of an additional ionization source 
which can produce strong [O III] and [N II] emission while producing very little Ha, H/3, 
[S II], and [O I], these large discrepancies effectively rule out cosmic rays as a strongly 
contributing ionization mechanism. However, we point out that, while they do not appear 
to be a sufficient ionization mechanism for the cool gas, cosmic rays may still be important 



in providing energy to the ICM, heating H2 (Ferland et al. 2009; Donahue et al. 2011), and 
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It is clear from Figures 14 16 that, individually, none of the models we listed at the 
beginning of this section can simultaneously explain all optical line ratios, their observed 
dependence on linewidths, and the total Ha luminosity. However, there is a very real pos- 
sibility that these filaments are, in fact, ionized by a combination of processes. Below, we 
discuss two such composite models which offer promising results. 



4-1.1. Composite Models - Star Formation in a Condensing ICM 
There is compelling evidence for both star formation and cooling fiows in the cool cores 



of many galaxy clusters. Modest cooling fiows are observed in the X-rays 
eraLl|2003l iVoigt & Fabianl 120041) and EUV (e.g.. 



(e.g., 


Peterson 


Bregman et al. 



2006), while the evidence for star formation is present at UV, optical, and IR wavelengths as 



we presented in §3.2. Thus, a composite model including self-ionizing, cooling gas and star 
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Fig. 16. — [N II] /Ha, [S II] /Ha, and [O I] /Ha line ratios as a function of the velocity 
dispersion for nuclei and filaments in our sample. Point types/colors are described in Figure 
3l In the left panels, the regions occupied by LIRGs and ULIRGs ( |Monreal-Ibero et al.|2010| 
are shown in grey and yellow, respectively. In the right panels, we show two-line fits to the 
data, with the Pearson R of each line in the upper left corner. Also shown in the upper 
left corner is the F-test probability, P2L; with which we can reject the null hypothesis of a 
single-line fit in favor of a two-line fit. 
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formation is both observationally motivated and straightforward to interpret. In order to test 
this model, we turn to Eq. [T| which predicts the Ha luminosity given an ICM cooling rate. 
We make the assumption that the cooling ICM surrounds the cool clouds and, thus, only 
half of the ionizing radiation is incident on the cool gas. We then compute the expected Ha 



luminosity on a cluster-by-cluster basis, using X-ray-derived cooling rates from McDonald 



et al. (2010) and McDonald et al. (2011a). Finally, we compare the predicted Ha luminosity 



to the measured values from the aforementioned references and determine the fraction of Ha 
luminosity in each cluster that is due to the condensing ICM. We further assume that this 
fraction is constant throughout the cluster and remove the contribution to each emission line 



flux from the cooling flow, assuming the predicted line ratios shown in Figure 14 
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Fig. 17. — Similar to figure [3} but with the expected contribution from the cooling fiow 
removed. By subtracting the expected fiux based on the X-ray cooling properties of each 
cluster, we find a better match between the optical emission line ratios in cool core cluster 
filaments and typical H II regions, with the exception of the [O I] /Ha ratio. 



In Figure [17] we again show various optical line ratios for the nuclei and filaments in our 
sample of 9 clusters, but now with the predicted contribution from the cooling fiow removed. 
We find that the [N II] /Ha and [S II] /Ha ratios are in better agreement with those for typical 
H II regions in most systems, with the exception of two systems classified as shock-heated by 



McDonald et al. (2011b), based on their Ha/FUV fiux ratios. However, the [N II]/Ha and 
[O I] /Ha ratios are still too high in most cases. Nevertheless, we feel that this composite 
model is a useful step forward since it introduces no additional free parameters and tries to 
account for a source of ionization which we know is present in cool core clusters. The optical 
line fiuxes due to the cooling ICM are tied directly to the X-ray-derived cooling rates and 
the observed Ha luminosities. 
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4-. 1.2. Composite Models - Star Formation in a Turbulent Cooling Flow 



Figure [7] shows that the hne-emitting gas in cool core clusters generally has a coherent 
velocity field. The combination of a bulk fiow on the order of ~ 200 km s~^ and line 
widths of a few times 100 km s~^ is reminiscent of the slow shocks combined with star- 
forming regions seen in large-scale winds ( Rich et aL||2010 ) and in luminous infrared galaxies 
(Rich et al. 2011). This composite model, proposed by Parage et al. ( 2010[ ) for the optical 
filaments in Centaurus A, produces slightly-elevated [N II] /Ha and [S II] /Ha, and highly- 
elevated [O I] /Ha, just as we observe in the cool cores of galaxy clusters. In Figure 18 
we show the predicted line ratios from a mix of star formation and radiative shocks from 
Allen et al. (2008). We find that a modest fractional contribution from shocks (~ 0-40%) 



in combination with star formation can produce optical line ratios that match the majority 
of our data. Further, if we account for self-ionization from the condensing ICM (§4.1.1), the 
required contribution from shocks is even lower (~ 0-20%). When combined with Figure [l6| 
this offers compelling support in favor of a composite model including both star formation 
and shocks. 
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Fig. 18. — Inspired by a similar plot in Rich et al. (2011), this figure shows the predicted 
emission line ratios for a mix of star formation and slow shocks over a range of ionization 
parameters and shock speed, assuming Z = Zq. The gray background points are data for 
filaments (stars) and nuclei (circles) in this sample. The upper panels show the measured 
fiux ratios, while the lower panels show fiux ratios with the contribution from a self-ionizing 
condensing ICM remove (§4.1.1, Figure 17). These models are able to reproduce the high 
[O I] /Ha ratios seen in the cool cores of galaxy clusters with a modest (~0-40%) amount of 
shocks. 



The observation of two populations in Figure 16, combined with the above confirmation 
that a mixture of shocks and star formation is sufficient to produce the observed line ratios. 
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suggests that we can divide our sample into star formation-dominated and shock-dominated 
subclasses. In Table 2, we summarize the properties of two subsamples: systems which have 
non-zero UV flux, low [N II] /Ha, and a weak or no correlation between the low ionization line 
ratios and velocity dispersion (star formation dominated; Abell 0478, Abell 0780, Abell 1795, 
Abell 2597, Sersic 159-03) and systems for which there is no detectable UV emission, high 
[N II] /Ha ratios, and a strong correlation between the low ionization line ratios and ve- 
locity dispersion (shock dominated; Abell 0496, Abell 1644, Abell 2052, NGC 4325). This 
table clearly demonstrates the differences between the star formation and shock dominated 
systems, in terms of their optical line ratios, optical line widths, UV flux and morphology, 
and nuclear line ratios. We caution, however, that the separation seen in Figure [16] is likely 
a selection effect and we would expect to see a continuum of systems with shock fractions 
varying from 0-100% in a larger, more complete sample. 

This scenario offers a straightforward explanation for, it seems, all of the observed 
properties of optical filaments in cool core clusters. The combination of shocks and star 
formation account for the soft X-ray, UV and excess blue emission, mid-IR emission, slightly 



elevated UV/Ha ratios (McDonald et al. 2011b), small linewidths, and elevated [O I]/Ha 



ratios. Including contributions from the cooling ICM reduces the need for shocks by roughly 
a factor of two, but is not necessary to match the data. Assuming that the H II regions 
are flowing along the filaments, they should be experiencing weak shocks due to interactions 
with the relatively stationary ICM. 



4.2. The Origin of Optical Filaments - Inflow or Outflow? 



The three most likely scenarios for the origin of the cool filaments in cluster cores are 
radial infall, buoyant outflow behind radio bubbles, and entrainment in radio jets. These 
models predict different kinematic signatures in the filaments which should be observable. In 
the classical picture of a rising bubble, cool material is entrained behind the bubble and rises 
to larger radius. Directly behind the bubble, the cool material should be moving away from 
the cluster center, while at small radius the gas which is trailing far behind the bubble cools 



and falls back into the bottom of the potential well (e.g., Churazov et al. 2001 Reynolds 



et al. 2005). This results in a stretching of the filament, with the inner and outer portions 



moving in different directions. The line-of-sight velocities of such a filament should show a 



change of sign along the filament, as is seen in the core of Perseus A (Fabian et al. 2003 



Hatch et al.||2006 ) . While we do not observe a change in direction along any of the extended 
filaments in our sample, this may simply mean that we are observing the filaments early 
on, at which time the bulk of the gas will be outflowing, or at a late stage, when all of 



- 32 - 



Systems: 



Star Formation Dominated 

A0478, A0780, A1795, A2597, 
S159-03 



Shock Dominated 

A0496, A1644, A2052, N4325 



Optical 

Emission-line 

Properties: 



< [N II] /Ha >™ 

< [S II] /Ha > nu 

< [O III]/H/3 >r 



: 0.99 ± 0.17 
0.89 ± 0.09 
= 0.25± 0.04 



< [N II]/Ha >fii = 0.85 ± 0.20 

< [S II]/Ha >fii = 0.61 ± 0.19 

[N II] /Ha oc 

[S II] /Ha oc ^0-74 



< [N II] /Ha 

< [S II] /Ha >„„c = 

< [O III]/H/3 > 

< [N II] /Ha = 

< [S II] /Ha = 

[N II] /Ha oc o-i-°6 
[S II] /Ha oc (71-32 



= 2.70 ± 1.24 
= 1.59 ± 0.54 
= 0.28 ± 0.12 

: 1.76 ± 0.33 
0.82 ± 0.21 



UV Properties": UV-Ha flux/morphology correlation Weak UV emission 



Nuclear Properties'*: Nuclear starburst 



LINER 



Shock Fraction*^ 



-0-40% 



-40-100% 



Table 2: Various properties of the two sub-samples which were defined based on how strongly 



the low ionization line ratios were correlated with the velocity dispersion (Figure 16). 



For systems which have HST far-UV imaging from McDonald & Veilleux (2009) and 



McDonald et al. (2011b). 



Based on optical line ratios in Figure^ 

Fraction of Ha luminosity produced by shocks, based on the composite star formation plus 



shocks models shown in Figure 18 



the cool gas is falling back on the BCG. Indeed, the extended filaments in Abell 2052 and 
Abell 2597 are coincident with the outer rims of the X-ray cavities, but have nearly-constant 
line-of-sight velocites. 

A second scenario for the extended, cool gas is that it is being uplifted by radio jets. 



This was recently proposed as the origin of the extended warm gas in Sersic 150-03 (Werner 



eraL][20lI| , but was also investigated much earlier for a sample of radio-selected galaxies by 



Baum & Heckman (1989). In this scenario, cool gas which has settled in the cluster core is 



entrained within radio jets during an episode of feedback. Removed from the direct infiuence 
of the AGN, this gas is allowed to cool further, resulting in the formation of young stars. 
The kinematics of the filaments in this case should be relatively smooth, with a constant 
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velocity direction dictated by the radio jet. This is consistent with our observations for the 
most extended systems, as discussed in §3.2. This process is hkely occurring in Sersic 159-03 



(Werner et al.||2011 


), and at small radii in Abell 1795 ( 


van Breugel et al. 


1984 


McDonald & 


Veilleux 


2009 


)• 





Gravitational freefall of a cool cloud may be the simplest of these three models, and is 
likely happening to some extent in most systems. Assuming a particle is in freefall through 
an atmosphere, the radial velocity should rise to a maximum speed at large radius where the 
atmosphere density is minimal, decelerate at intermediate radii as it experiences increasing 
drag, and then accelerate at small radii when gravity dominates once again. The particle 
should have the same velocity direction at all radii, unlike the buoyant bubble scenario. This 
characteristic shape is seen in the extended filaments of Abell 0496 and Abell 1795. Of the 
three aforementioned scenarios, this may be the most difficult to confirm, as there is no 
obvious secondary signature (i.e. radio emission. X-ray cavities). 

Without additional information, line-of-sight kinematics are insufficient to differentiate 
between the three aforementioned models. Seeking to identify whether optically-emitting 
gas is in the foreground or background of the BCG is critical to interpreting these velocities. 
Additionally, further work identifying the small-scale (arcsecond) orientation of radio jets in 
cluster cores and low surface brightness x-ray cavities would provide additional evidence in 
favor or against the two outflow models. 

5. Summary 

Utilizing a combination of narrow-band Ha imaging and carefully-positioned long-slit 
spectroscopy we have obtained optical spectra of Ha fllaments in the cool cores of 9 galaxy 
clusters. These spectra provide a wealth of information about the warm, ionized gas in 
cool core clusters, allowing us to provide the following new constraints on their origin and 
ionization mechanisms: 

• We flnd evidence for only small amounts of reddening in the extended, line-emitting 
fllaments, with the majority of fllaments having E{B — V) < 0.2. 

• We conflrm that the optical line ratios [N II] A6583/Ha, [O I] A6300/Ha, [S II] AA6716, 
6731/Ha are slightly elevated above the expectation for star-forming regions. We flnd 
that fllaments identifled as star-forming by their far-UV emission consistently have the 
lowest [N II] /Ha ratios. 

• For the complete sample there is a weak correlation between the low-ionization line 
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ratios and line widths. This correlation is stronger for systems which appear to have 
a larger contribution from shocks than from star formation, based on their optical line 
ratios and far-UV/Ha flux ratios. 

The optical line ratios are able to rule out collisional ionization by cosmic rays, ther- 
mal conduction, and photoionization by both the ICM X-rays and AGN as strong 
contributors to the ionization of the warm gas in both nuclei and filaments. 

A composite model of a self-ionizing cooling plasma and star formation is able to 
better reproduce these line ratios than star formation alone, but still under-predicts 
the low-ionization ratios in most systems. 

The data are fully consistent with a combination of slow shocks and star formation, as 
seen in superwinds and merger remnants. This scenario explains the kinematics, line 
ratios, and multi-wavelength properties of the filaments and nuclei. We find typical 
shock fractions of 0-40% in systems with signatures of star formation (i.e. far-UV 
emission), and 40-100% in systems with no observational evidence for star formation. 
Including self-ionization from the cooling plasma to this model reduces the necessary 
contribution from shocks by roughly a factor of 2 while still matching the data. 

The nuclei of Abell 0496, Abell 1644, and Abell 2052 are identified as LINERs. The 
remaining 6 nuclei in the sample share similar properties to the star-forming filaments 
and likely have very little contribution from an AGN. 

In the central (r < 10 kpc) regions of several cool core clusters, we see kinematic 



signatures of rotating. Ha structures, consistent with earlier work (e.g., Baum et al. 



1992). 



The line-of-sight velocities of the most extended filaments peak at ~300 km s~^. We 
do not observe a change in direction along any filaments. We provide evidence drawn 
from the literature for entrained gas in radio jets (Sersic 159-03, Abell 1795), cool 
gas uplifted by buoyant bubbles (Abell 2052, Abell 2597), and gravitational freefall 
(Abell 0496, Abell 1795). 

The line profiles of the ionized filaments have widths on the order of FWHM ~ 200 
km s~^. This is consistent with velocity broadening due to slow shocks. In general, 
filaments that are thin and extended have narrower emission lines than those that 
appear disturbed. 

The line widths in the filaments decrease with distance from the BCG center. If the 
cool filaments have properties which reflect the ICM environment, this provides an 



- 35 - 



estimate of the radial velocity dispersion profile in the ICM. These velocities widths 
are consistent with those measured for the ICM using X-ray grating spectrometers 



e.g., Sanders et al. 2011) and are much smaller than the stellar velocity dispersions 



in the BCG. 



These new data, coupled with our previous work ( [McDonald etal^lMol [Mla||bl ) , sup- 
port a scenario where some fraction of the cooling ICM in the cluster core has been allowed 
to collapse and form stars. While in bulk motion (either infall or outflow), the cool gas is 
experiencing shocks, producing the observed line widths of ~ 200 km s~^, which are signif- 
icantly higher than those measured in typical H II regions. Further, the cooling ICM may 
contribute a non-negligible fraction of self-ionization as it cools through EUV and UV transi- 
tions. In filaments which appear to be star-forming based on their FUV and IR emission, the 
combination of shocks and self-ionization yield 0-40% of the Ho; luminosity in the filaments, 
with photoionization by young stars providing the additional 60-100%. Systems without 
any accompanying FUV emission appear to be fully consistent with ionization by radiative 
shocks. While satisfactory, this scenario is necessarily incomplete. A complete modeling of 
the multiphase plasma using all available data from X-ray through radio, rather than focus- 
ing on a single energy regime, should allow for a more robust solution to this long-standing 
problem, while combining the warm gas kinematics with deep X-ray and radio studies may 
shed new light on the motions of this gas. 
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Appendix B 

Below, we provide examples of tlie quality of the spectra used in this paper. A total of 
four spectra are shown, two from the shock-like emission in NGC 4325, and two from the 
star-forming filaments in Abell 1795. These two systems cover a broad range of signal-to- 
noise, from very high in the centers to our detection limits at large radius. Each spectrum is 
divided into 5 segments in order to clearly show the quality of the fits to the various optical 
emission lines. Each row represents a different spatial bin along the slit. Vertical dotted 
lines represent the systemic velocity of the BCG. 
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Fig. 19. — NGC 4325: The broad emission lines and high [N II] /Ha ratios suggest that 
shocks may be contributing a significant amount of the total ionization in the warm gas. 
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Fig. 19.— 



Continued (Slit #2). 
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Fig. 20. — Abell 1795: The narrow emission lines and small [N II] /Ha suggests that emission 
in these filaments is mostly due to star formation. The high [O I] luminosity requires a small 
fractional contribution (~20%) from shocks. 
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g. 20.— Continued (Slit #2). 



